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Abstract. Recently, 59 low-luminosity object transits were reported from the Optical Gravitational Lensing Experiment 
(OGLE). Our follow-up low-resolution spectroscopy of 16 candidates provided two objects, OGLE-TR-3 and OGLE-TR-10, 
which have companions with radii compatible with those of gas-giant planets. Further high-resolution spectroscopy revealed 
a very low velocity variation (< 500 m s _1 ) of the host star OGLE-TR-3 which may be caused by its unseen companion. An 
analysis of the radial velocity and light curve results in M < 2.5 Mj up , R < 1.6 Rj up , and an orbital separation of about 5 Rq, 
which makes it the planet with the shortest period known. This allows to identify the low-luminosity companion of OGLE-TR-3 
as a possible new gas-giant planet. If confirmed, this makes OGLE-TR-3 together with OGLE-TR-56 the first extrasolar planets 
detected via their transit light curves. 

Key words, binaries: eclipsing - stars: individual: OGLE-TR-3 - stars: low-mass, brown dwarfs - stars: planetary systems 



1. Introduction 

The detection of planets outside the solar system was a long- 
standi ng goal of astronomy. Af t er the first successful d etec- 
tions dWolszczan & Fraill 1 19921 iMayor & Ouelozl Il995l) . the 
search with various methods ("see ISchneiderf l200ll for an 
overview) was largely intensified. Out of the 104 presently 
known planets, 100 have been detected with Doppler- velocity 
measurements of the planets' host stars. All these plan- 
ets were found around solar-like stars. The other four 
are planets around pulsars and were found by periodic 
pulse-modulation measurements (see The Extrasolar Planets 
Encyclopaedia http://www.obspm.fr/encycFcatalog.html for an 
up-to-date overview). 

Until very recently, no planet had yet been found by pho- 
tometric monitoring. HD 209458 was the only system which 
has an orbital inclination that allows th e measurement of 
the eclipse of the host star by the planet JCharbonneau et alJ 
2000; lHenrv et all2 000). This planetary compani on was, how- 
ever, known before from Doppler measurements (M azeh et alJ 
2000f). The current paucity of extrasolar planets found via the 
transit method is deplorable for several reasons. The transiting 
systems provide more reliable parameters. With the known in- 
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clination, the radial velocity variation provides a precise mass 
determination and not only a lower mass limit. Together with 
the radius, the density of the planet can additionally be de- 
rived. The transit also opens, at least in principle, the pos- 
sibil ity to explore the pla net's atmosphere during the eclipse 
dCharbonneau et alll2002l) . Furthermore, transits are the most 
promising way to detect earth-like planets around solar-like 
stars and it is therefore the method of choice for several ground 
and space based projects. 

Recently, 59 transiting planet candidates were announced 
by the O ptical Gravitational Len sing Experiment (OGLE) con- 
sortium ( Uda lski et alJ l2002albl) . These candidates were ex- 
tracted from a sample of about 5 million stars observed during a 
32-day photometric monitoring. In a sub-sample of 52000 stars 
with a photometric accuracy better than 1.5%, these 59 candi- 
dates exhibit light curves indicating the presence of a transit- 
ing low-luminosity companion. From t his sample we obtain ed 
low-resolution spectra for 16 objects ( Dr eizler et alJl2002l) in 
order to derive spectroscopic radii of the primary stars and to- 
gether with the eclipse light curve also the companion radii. 
While 14 low-luminosity companions could be identified as 
M stars, two objects, namely OGLE-TR-3 and OGLE-TR-10, 
have companions with radii compatible with those of gas-giant 
planets. OGLE-TR-10 has been further investigated by Hatzes 
et al. (ESO Programme 268.C-5772). We will present in this 
paper our VLT follow-up observations for a dynamical mass 
determination of the companion of OGLE-TR-3. 
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Very recently, iKonacki et alJ J2003t) claimed the verifica- 
tion of the planetary nature of OGLE-TR-56. They also investi- 
gated OGLE-TR-3, OGLE-TR-10, OGLE-TR-33, and OGLE- 
TR-58. While the interpretation as planetary host star for none 
of these objects is conclusive from their work yet, they regard 
OGLE-TR-3 as the result of blending and grazing eclipse in 
contrast to our interpretation. 

In the following, we describe our observations and data re- 
duction (Sect. 2). The analysis is presented in Sect. 3 and the 
results are discussed in Sect. 4. 
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Fig. 1. From top to bottom: A small section of the co-added 
OGLE-TR-3 spectrum and the individual spectra of the ten ob- 
serving blocks (1 h observing time each) ordered by JD as listed 
in Tabled 



2. Observations and Data Reduction 

The aim of the observations is a detection of the reflex motion 
of OGLE-TR-3 due to the low-lu minosity companion caus- 
ing the observed eclipse light curve (Udalski et al. 2002a). The 
short orbital period (P= 1.1899 days) in combination with a 
spectroscopic mass of 1 . 1 M© (breizler et alJ l2002h . a low- 
mass star (M > 0.08 Mq) would produce a radial velocity 
variation of the order of 10km s _1 . A discrimination between 
a stellar and sub-stellar object therefore needs only moderate 
instrument requirements. We observed in Service Mode with 
UVES mounted at the VLT Kueyen. We used the red arm of 
the spectrograph with a central wavelength of 5800 A, a slit 
width of CK'7, and the CCD detectors without binning. Image 
sheer #1 was inserted in order to obtain the full spectral resolu- 
tion of 60000 independently of the seeing. 

We were granted one night of Directors Discretionary Time 
which we spread in ten observing blocks of each 3x20min 
exposures to minimize the velocity smearing due to the earth's 
motion. We asked for five consecutive nights with each two 
observing blocks around the meridian passage of OGLE-TR-3 
in order to sample the orbital period equidistantly. Due to bad 
weather and other targets of opportunity this schedule could 



not be kept and observations were finally spread out over about 
one month (for details, see Table[0. 

The spectra were extracted using procedures of the UVES 
context of the ESO-MIDAS software package. The routines 
provide correction for bias, flat-field, and background level. 
Wavelength calibration was performed with ThAr spectra ob- 
tained immediately after each observing block. According to 
the instrument handbook this provides a precision to about 
100ms- 1 . 

From the exposure time calculator we expected a maximum 
S/N of about 15 in a 20min exposure. The achieved S/N was 
significantly lower (around 10). We therefore had to co-add the 
three individual frames of one observing block to each one sin- 
gle exposure. As an example, we show a small section of the 
co-added spectrum and of all spectra from ten observing blocks 
(Fig.0 overlay ed with a synthetic spectrum from our spectral 
analysis (see below). 

3. Analysis 

3.1. Spectral Analysis 

An inspection of the VLT-UVES spectra of OGLE-TR-3 re- 
veals three spectral components. The earth's atmosphere com- 
ponent, consisting of the night sky emission and telluric lines, 
are easy to separate since they do not follow the radial velocity 
shifts due to the earth's orbital motion. Besides the stellar ab- 
sorption lines, there is a third component, mainly visible in the 
NaD doublet at about 5890 A, which has a constant velocity 
offset of 20km s _1 . The second absorption component could 
be caused by a stellar companion of OGLE-TR-3 or by inter- 
stellar absorption. If this were due to a binary companion we 
could use the difference in radial velocity and interpret it as or- 
bital velocity. With the parameters of OGLE-TR-3 this would 
result in an orbital period of about 4 years. This component 
can therefore not be responsible for the observed 1.1 899 day 
eclipse period. If identified as binary companion we would ex- 
pect to see other spectral features at that velocity shift, if it is 
also a late F to early G star as OGLE-TR-3. This is, however, 
not the case. The possible stellar companion could alternatively 
be of later spectral type, down to late K. It would then be about 
a factor of two fainter. Since NaD is among the strongest spec- 
tral features for those spectral types, this would explain why we 
cannot detect other lines. We should, however, detect the MgH 
band around 5100 A, which is not the case. 

As demonstrated in Fig.|2] we can fit the NaD line nicely 
with an interstellar absorption caused by a column density of 
5 • 10 15 cm~ 2 neutral Na atoms. This is a very plausible value 
since OGLE-TR-3 lies about 3 kpc in the direction of the galac- 
tic bulge. 

We also investigated the possibility of blending, first due to 
a second, equally luminous star at the velocity of OGLE-TR-3. 
This hypothetical third object would then be a late F to late G 
star, since an earlier spectral type would dominate the spectrum 
and a later one would hardly contribute to the total brightness. 
The eclipse could then result from a late K or early M star or- 
biting one of the G stars, the shallow eclipse would just be a 
fake due to the blend (G+[G+M] scenario in the following), 
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Fig. 2. Three normalized VLT-UVES spectra of OGLE-TR-3, obtained at JD-2452400= 62.61590, 86.57832, and 89.56877 (from 
bottom to top), around the NaD doublet. The fit consists of three components: A synthetic spectrum shifted by the measured radial 
velocity (blue wing), a Gaussian emission component from the earth's atmosphere, and an interstellar absorption of a column 
density of 5 • 10 15 cm~ 2 neutral Na atoms for the red wing. The velocity of the interstellar component has a fixed offset of 
+20kms -1 relative to the stellar component. 



as suggested bv lKonacki etal] §003). We therefore carefully 
inspected our spectra in order to detect asymmetries in the line 
profiles. As demonstrated in Fig.[^with the overlayed synthetic 
spectrum, all individual spectra appear symmetric and the co- 
added spectrum does not show any signs of additional broad- 
ening. The effect, if present, is expected to be directly visible 
since an M star companion to a G star in the short orbital pe- 
riod of OGLE-TR-3 would result in a radial velocity variation 
of the order of 10km s _1 corresponding to a one pixel shift on 
the detector. This would also be visible in the co-added spec- 
trum as additional broadening of the same order of magnitude. 
This is not the case (see below). For a more quantitative inves- 
tigation, we fitted 17 strong and unblended spectral lines with 
Gauss profiles and calculated the difference between the obser- 
vations and the Gauss fit for the red and blue wings separately. 
The total difference is shown in Fig. [5] Any asymmetry would 
result in systematic offsets between the red and blue line wings. 
Additionally we cannot find a trend with the orbital period. Our 
current data set provides no hint for blendi ng and we th erefore 
regard the scenario proposed by iKonacki et alJ J2003I) as un- 
likely in case of OGLE-TR-3. 

The second blending scenario is a close M+M star binary 
orbiting at 2x1. 1899 days blended by a G star with a large 
enough distance towards the M star binary so that it has an un- 



detectable orbital velocity (G+[M+M] scenario in the follow- 
ing). We inspected our UVES spectra in order to detect spectral 
features from a possible M star. There is no evidence for it but 
since the contribution of two M stars is only of the order of 3% 
in the red, we cannot rule out this completely. The light curve 
from such a system is discussed below (Sect. 3.3). 

The third blending scenario is a line of sight blending: A 
GO star is contributing most of the light while an eclipsing bi- 
nary system is situated at a much larger distance (G+[X+Y] in 
the following). Its contribution to the spectrum must be (very 
conservatively) less than about 30% because we would detect 
it otherwise, as discussed above. Its influence on the light curve 
is also discussed below (Sect. 3.3). 

We fit the co-added spectrum of OGLE-TR-3 by compar- 
ing it to a PHOENIX (version 13.00) spherical model atmo- 
sphere grid and synthetic spectra calculated at the resolution 
of the observed spectrum (Figs.^ and 0}. The model atmo- 
sphere grid covers a wide range of effective temperatures, grav- 
ities, and abundances . It includes the lates t updates with re- 
spect to the NextGen jHauschildt et aljfl 999) and Dusty /Cond 
lAllard et al J 1200 ll) grids. For simplicity, we have used LTE 
models, since NLTE effects in the relevant parameter range can 
be expected to be small and were thus neglected for the initial 
fits. The updates include improvements in the equation of state 
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Fig. 4. Four 50 A sections of the co-added spectrum of OGLE-TR-3 compared to a theoretical model: T g = 6100 K, log g = 4.5 
(cgs), and sub-solar abundances ([M/H] = —0.5). 



data used and better molecular line lists, however, the latter are 
not important for OGLE-TR-3. 

The synthetic spectra grid covers a range of 6000 K< 
T e H < 6200 K and 4.0 < log g < 5.0 for solar abundances and 
[M/H] = —0.5, in agre ement with our earli er spectral classifi- 
cation as F9/G0 V star dDreizler et alJl2002h . Synthetic spectra 
were computed with 0.01 A steps and rotationally broadened 
for u rot sini = 2kms~ 1 . Again, this is a clear indication that 
no third object as in the G+[G+M] scenario is present, as dis- 
cussed above. In addition, the spectral resolution was down- 
graded to 0.02 A before comparison with the data. The ob- 
served spectra are not flux calibrated, so the synthetic spec- 
tra were normalized to the continuum. With this procedure, we 
find that the best fit parameters for OGLE-TR-3 are T e ff = 
6100 ± 200 K, log g = 4.5 ± 0.5 and abundances between so- 
lar and [M/H] = —0.5. We did not attempt to improve the fit 
by changing the abundance pattern, but inspection of the com- 
parisons shows that this could lead to improved fits. Models 
outside the parameter range do not lead to improved fits, and 
were not considered further. In Fig.|5)we compare these param- 
eters with stellar evolution models I Schaerer et al Jl 19931) and 
derive the stellar mass and radius. The spectral analysis cor- 
roborated the stellar parameters of the primary derived earlier 



(R = 1.05 R©, M = 1.1 M B iDreizler et aljEool . However, 
our fit favors slightly lower parameters, R = l.O^o 2 R© an d 
M = l.O±g;?M . We adopt R= 1.0 R Q and M — 1.0M Q in 
the following radial velocity and light curve analysis. 



3.2. Radial velocity variation 

For the determination of the radial velocity variation of OGLE- 
TR-3 as reflex motion due to the unseen companion, we 
used the first spectrum as reference and determined the cross- 
correlation in velocity space relative to the other nine. We iden- 
tified twelve 50 A intervals within the available spectral range 
(4800-6800A) with little contamination from telluric lines. 
Additionally, we had to avoid the spectral region between 5750 
and 5850 A where the spectrum falls onto the gap between the 
two CCDs. The cross-correlation function is approximated by 
a Gauss profile to determine the velocity shift at sub-pixel pre- 
cision (see Fig. [6}. The median of the twelve segments is then 
taken as the velocity shift of the spectrum. The standard de- 
viation is used as quality attribute of the measurement and di- 
rectly reflects the S/N of the spectra. Afterwards we co-added 
all spectra in velocity space and repeated the procedure, now 
with the co-added spectrum as reference. Reassuringly, neither 
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Fig. 5. The parameters of OGLE-TR-3 compared to stellar evolution models JSchaerer et alJl993l) . The evolution tracks as well 
as the lines of constant radii are labeled in solar units. 
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wings of observed spectral lines and Gauss profiles with 1 a 
error bars. The spectra from JD=86. 53247 and 86.57832 are 
omitted due to very large uncertainties. No indication of line 
asymmetry is present in our data set. 



Table 1. Observing dates, orbital phase of OGLE-TR-3, mea- 
sured radial velocity relative to a co-added spectrum, earth's 
orbital and earth's rotational velocity components towards 
OGLE-TR-3, radial velocity of OGLE-TR-3, and standard de- 
viation of the measurements. The radial velocity column pro- 
vides the measured velocity corrected for earth's motions mi- 
nus the median of the 10 values. 



JD-2452400 


phase 




radial velocity [km s x ] 








spec. 


orbit 


rot. 


TR-3 




62.61590 


0.1754 


0.06 


-8.32 


0.11 


0.151 


0.12 


62.66215 


0.2143 


0.15 


-8.34 


0.01 


0.116 


0.17 


63.63365 


0.0307 


0.51 


-8.80 


0.07 


0.072 


0.10 


63.67882 


0.0687 


0.56 


-8.82 


-0.03 


-0.003 


0.16 


66.61967 


0.5400 


2.11 


-10.20 


0.08 


0.290 


0.25 


66.66758 


0.5804 


1.95 


-10.22 


-0.02 


-0.004 


0.13 


86.53247 


0.2741 


10.11 


-18.65 


0.15 


-0.094 


0.58 


86.57832 


0.3127 


10.32 


-18.67 


0.05 


0.000 


0.64 


89.52445 


0.7885 


11.14 


-19.77 


0.15 


-0.192 


0.44 


89.56877 


0.8257 


11.13 


-19.79 


0.06 


-0.306 


0.55 



the cross-correlation with the co-added spectrum nor with a 
model spectrum does change the results significantly. 

Using the period (P=1.1899day s) and epoch of the mid 
eclipse (JD=2452060.22765) from lUdalski et ail d2002al) . we 



calculated the phasing of each mid of exposure . For each ob- 
servation we also calculated the velocity components of the 
earth towards OGLE-TR-3 using the BARYVEL routine avail- 
able in the Interactive Data Language (IDL) astrolib from 
Goddard Space Flight Center. It should be noted, that this rou- 

1 The slightly different values given at the OGLE-web site 
http://sirius.astrouw.edu.pl/~ogle result in the same conclusions. 
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Fig. 6. An example of the cross-correlation. Co-added spec- 
trum versus the one from JD=2452486. 53247 in the interval 
[5350 A- 5400 A] (++++) and a Gauss fit (solid line). 
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Fig. 8. Measurements of the velocity shift of the night sky 
emission lines in the spectra of OGLE-TR-3. Error bars indi- 
cate 1 a deviations. The solid line is the best fit from our \ 2 - 
minimization. 



tine provides the velocity in a right-handed coordinate system 
with the +X axis toward the vernal equinox, and +Z axis to- 
ward the celestial pole. A positive velocity component towards 
an object therefore produces a blue shift whereas the spectro- 
scopic definition for the sign of the velocity is the opposite. We 
also calculated the velocity components due to the earth's ro- 
tation, which is, however, of minor importance. In Table[Jwe 
summarize our results for each exposure. 

An independent period determination is impossible from 
our few scattered measurements. Hence we mapped the radial 
velocity curve t o the orbital phase (Fig .|7t using the period and 
zero point from Udals ki et alJ (|2£)02a). These data were fitted 
with a sine curve with an offset, the velocity amplitude, and a 
phase shift as free parameters. We can use the fitted phase as 
an indication that the period is present in our data. As expected, 
the derived radial velocity is indeed zero at mid eclipse. 
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Fig. 7. Measurements of the radial velocity of OGLE-TR-3. 
Error bars indicate 1 a deviations. The solid line is the best fit 
from our ^-minimization, the dashed lines are the 3 a devia- 
tions of the amplitude. 



First, we used an inversion based x 2 -minimization routine. 
Fig. [7J shows the measured radial velocities with error bars of 
1 a overlayed with the best fit (solid). We also show the veloc- 
ity curve for 3 a variation (dashed) in the amplitude (see also 
Table|2}- The resulting amplitude of about 100ms -1 is smaller 
than the median 1 a error, we therefore investigate the signifi- 
cance of the result in more detail. 



Table 2. Results of the radial velocity curve fits. Left block: 
Results from the inversion based x 2 -minimization. Right 
block: Results from the genetic algorithm optimization. The 
blocks give the best fit value and lower and higher limits. Note 
that the 3 a deviation of the pikaia fits are smaller than 1%. 



parameter 




x 2 




pikaia 




best 


-3 a 


+3 a 


best 


amplitude [kms _i ] 


0.10 


0.08 


0.13 


0.12 


phase 


0.02 


0.99 


0.04 


0.94 


offset [km s~ ] 


0.05 


0.04 


0.06 


0.05 



We used night-sky emission lines which do not show veloc- 
ity variation due to the earth's motion. The resulting "velocity" 
variation (Fig. [8} reveals a fit with vanishing amplitude, as ex- 
pected. It also demonstrates that we recover the velocity preci- 
sion of UVES (« 100 m s _1 ). This test is an indication that the 
derived velocity amplitude is indeed due to intrinsic variations 
rather than to measurement errors. 

Additional ly, we utilized the genetic algorithm pikaia 
dCharbonneaul !l995) to confirm the three fit parameters. The 
genetic algorithm starts from randomly distributed parameter 
sets and increases the quality of the fit according to selection 
criteria adopted from the biological process of evolution. It 
therefore provides the global x 2 -minimum much more reliably 
than classical x 2 -minimization. In order to obtain a measure- 
ment for the reliability of the result, we ran pikaia with 1000 
different initializations of the starting population. The deviation 



Dreizler et al.: OGLE-TR-3: A Possible New Transiting Planet 



1 .0 =~ 
0.8 Er 
'.6 EE 

■A — 



.1 

| 

1 

> 



.3e-5 



0.1 170 



= 

8 E- 

6 Er 

4 E- 

2 E- 
o EE_ 



0.1 1 75 



0.1 u 



RV observation 



0.1 185 
[km/s] 



1" 



0.1 190 



0.1 195 



•is 

_o 

1 

•" 

CD 

en 

o 

a 



0.9380 



0.9385 



0.9390 
phase observation 



3.4e-5 



0.9395 



0.9400 





-0 

-0 



2 — 
3- 
2 — 



0.0515 



0.0520 



0.0525 
offset observation [km/s] 



5.6e-5 



0.0530 



0.0535 



Fig. 9. Fit results of synthetic white noise curves versus fit results of the OGLE-TR-3 radial velocity (RV) curve 
While the synthetic data (ordinate) spread over a wide interval, the real data (abscissa) are confined to a narrow 
ratio between the distribution widths of the fits is indicated in the plots. 



using pikaia. 
range only. The 



of the final fit parameters is very small (Table|2ji and they are 
consistent with the one of the classical \ 2 -minimization. 

As a third test, we synthesized 1000 white noise curves 
with the variance of that from the radial velocity curve and with 
identical timing and fitting weights. These data sets were fitted 
with pikaia. In Fig.|5]we compare the results of the pikaia 
fit of the radial velocity curve with those from the simulations. 
While the parameters of the radial velocity variations are recov- 
ered with a very small scatter, the noise simulations produce 
fit results over a large parameter space. While the fitted ampli- 
tudes and offsets are typically of the order of the mean standard 
deviation of our measurements (« 200 ms -1 ), the phase fits 
are equally distributed over the interval [0, 1 ] . As measurement 
for the reliability of the fit, we provide the widths of the param- 
eter interval for the radial velocity fits divided by the width of 
the simulation fits. 

Concluding, we can fit the observed radial velocity curve 
with an amplitude of about 120 m s _1 . We provided several in- 
dications that this fit indeed represents a true reflex motion of 
the late F/ early G type planetary host star due to the unseen 
companion. This contradicts the G+[M+M] scenario, where the 
G star would have a very small radial velocity variation at a 
much larger orbital period. 

3.3. The light curve revisited 

For a consiste ncy check, we fitted the light curve of OGLE-TR- 
3 published by Udalski et al. ( 2002a) with our derived parame- 
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Fig. 10. Light curve of OGLE-TR-3 from llJdalski et ail 
( 2002a): Top and middle: I band observed (++++) and synthetic 
light curve (solid line). Bottom: Predicted H band light curve at 
a smaller scale. The primary eclipse is shown in the top panel, 
the secondary eclipse in the lower two. 

ters (R — 1.0 R©, M = 1.0 M Q ). If we estimate the tempera- 
ture of the secondary due to irradiation, we derive about 2 1 00 K 
for a 6100 K primary and 5 R© separation. We utilized the bi- 
nary eclipse simulation program nightfall (R. Wichmann, 
Landessternware Heidelberg, Germany) which calculates syn- 
thetic light curves taking into account the distortion of the stars 
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in Roche geometry. As can be seen in Fig. [TO] we can nicely 
reproduce the observed I-band light curve with an inclination 
of 90°, a primary radius of R = 1.25 R Q and a secondary 
radius of 0.14R Q (80% filling factor of the Roche potential 
of a 0.0006 M Q planet). The radial velocity of the primary is 
131 kms -1 for such a system, fully consistent with our deter- 
mined value. The rotational velocity of the primary is 2 km s _1 
which is far below synchronization («5%). Due to the scat- 
ter in the photometry a simultaneous fit of the inclination and 
the secondary radius is not very well constrained. An inclina- 
tion of 85° requires a secondary radius of 0. 16 Kq, 80° require 
0.17Rq. The formal fit quality, however, favors an inclination 
very close to 90°. Wh ile we derive a radius ratio consistent 
with that from Udalsk iet alJ (|2002a), our absolute radii are 
slightly smaller than their values (i? x =1.48 R Q , i? 2 =0.18 R Q ), 
but within our uncertainty. A very small secondary eclipse 
might be present in the current photometry. Whether this is real 
remains to be checked with much more precise photometry. 
Fig.^|also shows the predicted H-band light curve of the sec- 
ondary eclipse which is of the order of 1 mmag. It should there- 
fore be possible to detect the planet directly with high precision 
IR photometry. We also investigated the effect of the distortion 
of the primary due to the gravitational field of the secondary. 
This effect becomes visible in a sin(2P)-variation outside the 
eclipses. We fitted both our model solution as well as the ob- 
servation and find a 0. 15 mmag and 0.3±0.3 mmag variation (3 
a error), respectively. Our solution is therefore consistent with 
the observations. 




0.3 0.4 0.5 0.6 0.7 

phase 



Fig. 11. Light curve of OGLE-TR-3 from lUdalski et all 
<2002~1: I band observed (++++) and synthetic light 
curve (solid line). The photometric data are folded with 
2xl.l899days in order to simulate the G+[M+M] scenario. 
The eclipsing system (primary minimum top, secondary min- 
imum bottom panel) consists of two M2 stars; the third light 
contribution is that from a G star. 

While our radial velocity measurements are in favor of the 
interpretation of OGLE-TR-3 B as planet, the light curve might 
be explained in the G+[M+M] scenario. The eclipse of the two 
M2 stars diluted by the GO star contribution would result in 
a shallow eclipse of the observed magnitude difference. It re- 
quires that the two M stars are nearly identical. The measured 
period would correspond to half an orbital phase and the true 
period of the system would then be 2 x 1.1 899 days. A period 
search in the original data, however, results in the largest peak 



at the period of 1.1 899 days. If we simulate the light curve of 
this system including the third light contribution, we obtain a 
reasonable fit (Fig. ll 1> . which is, however, not as good as in the 
case of a transiting planet (Fig. llOt . The eclipses are too nar- 
row, the formal fit error is therefore larger than for the planet 
transit light curve. An increase of the stellar radii to increase 
the eclipse widths would also influence the depths since the bi- 
nary would contribute more to the total light. We therefore re- 
gard this interpretation for the OGLE-TR-3 system as unlikely. 
The G+[M+M] scenario can easily be checked, since the M star 
contribution to the total light is strongly wavelength dependent. 
Hence, we would expect a strong color effect in the light curve. 

The G+[X+Y] blending scenario would also require a 
false detection of the radial-velocity variation. Furthermore, 
we checked possible blends for their consistency with the light 
curve. Assuming that the secondary is much fainter than the 
primary (no detectable secondary eclipse), we can pick a spec- 
tral type of the secondary with the corresponding radius and 
mass for a main sequence star. Unblended light curves for vari- 
ous contributions from a less distant GO star are reconstructed, 
the depth of the primary eclipse then provides the radius of 
the primary. The corresponding mass together with the period 
is used to derive the orbital separation and velocity (Kepler's 
third law, assuming circular orbits). The time between first and 
fourth contact can then be obtained from the orbital velocity 
and the sum of the radii, the time between second and third 
contact from the orbital velocity and the difference of the radii. 
In Fig. [21 these eclipse times are plotted in units of the cor- 
responding times obtained from the observed light curve. A 
consistent solution is only possible in a very limited parameter 
range: The foreground GO star provides about 80% of the total 
light, the eclipsing binary system is then a F9 secondary and 
a B9 primary. However, this scenario can be ruled out when a 
theoretical light curve is calculated with these parameters. The 
close vicinity of the two stars would result in an observable 
sinusoidal variation. Additionally, the light curve can only be 
reproduced very roughly with this solution. We therefore reject 
the G+[X+Y] scenario based on our current observational data. 



4. Summary and Discussion 

Summarizing, our data favor the interpretation of OGLE-TR-3 
as planetary host star and we have detected the radial veloc- 
ity variation due to the reflex motion of the unseen planetary 
companion. 

From Kepler's third law together with momentum conser- 
vation the compani on mass can be determined from the period 
(1.1 899 days, Udal ski et alJl2002al) . the mass of the primary 
(1.0 M Q , this paper), the velocity amplitude (120ms _1 , this 
paper), and the inclination, which must be close to 90° (see 
above, in the following we will use 90°). 

The upper limit of the derived radius (R = 0.17 Rq) of the 
unseen companion would be consistent with an M star down 
to a gas-giant planet. With our VLT observations, however, 
we can clearly rule out an M star since a 0.08 M Q star would 
cause a 15kms _1 reflex velocity variation. We can also rule 
out a brown dwarf companion since a brown dwarf at the lower 
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Fig. 12. Predicted eclipse durations for blending scenarios as 
function of the spectral type of the secondary (first to fourth 
contact: full line; second to third contact dashed line) in units 
of eclipse durations of OGLE-TR-3. Contribution from a fore- 
ground star with 70, 80, 90, 95, 98% is labeled with crosses, 
stars, squares, diamonds, and triangles respectively. 



limit for deuterium burning (0.013 M Q ) would still result in 
2.5 kms" 1 velocity amplitude. As a conservative upper limit 
for velocity variations (3 a) we derive 500 m s" 1 which trans- 
lates into 0.0026 M . If our fit is regarded as reliable, we obtain 
a mass of 0.0005 to 0.0006 M or about half a Jupiter mass. 
Due to the small period, the separation of the two objects is 
only about 5 R which is however still a factor of two above 
the Roche stability limit. The companion's mean mass den- 
sity is comparable to the one of HD 209458 B. If we use our 
0.0005 M Q and 0.014 R Q , it amounts to 0.25 gem" 3 , which is 
about 20% of the mean solar density. 

The shortest known orbital period of all extrasolar planets 
and the very close orbit of OGLE-TR-3 raises the question of 
the influence of tidal effects between the star and planet. As the 
star appears to be rotating much slower than the planetary or- 
bital period, any tides raised by the planet on the star will lead 
to an additional shrinki ng of the o rbit. Following the argmu- 
ments of iGoldreich & Soteil i 19661) we estimate the timescale 
a I a for a change in the semi-major axis of the planet, using 
the present orbital parameter. Using a value of 10 7 for the dis- 
sipation factor Q* of the star we obtain a damping timescale 
of about 2 ■ 10 9 yrs. As the Q value for stars is quite uncertain 
and may easily be higher, we may conclude that a planet is just 
about able to survive in such a close orbit for the lifetime of the 
star. 

We therefore conclude that we confirmed the low- 
luminosity companion of OGLE-TR-3 as an extraso- 
lar planet (Mi=l.OM , Af 2 =O.OOO5M and i? 2 =O.14R ) 
with quite similar properties compared to HD 209458 
(Mi=1.05M Q , M 2 =O.OOO66M , i? 2 =O.147R ) and OGLE- 
TR-56 (Mi=l.O4M ,M 2 =O.OOO86M Q ,i? 2 =O.13R ) regard- 
ing the parameters of the primary star and planetary compan- 
ion. The system of OGLE-TR-3 is slightly more extreme re- 
garding the separation. 



Finally, if confirmed, our results for OGLE-TR-3, together 
with those for OGLE-TR-56, make these objects the first two 
extrasolar planets detected via the transiting method. A further 
confirmation of the planetary nature would benefit from sev- 
eral supporting observations: Optical B,V,R photometry to ex- 
clude a color effect expected in the blending scenarios, IR spec- 
troscopy to exclude an M star contribution, H band photometry 
to detect the predicted very shallow secondary eclipse, and fi- 
nally, more radial velocity measurements to detect the orbital 
period in the radial velocity data directly. 
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